Nuclear effects in high-p^ production of direct photons and neutral mesons 



in 
o 
o 

(N 

G 
3 



> 

m 
o 
o 
\o 
o 
in 
o 

">< 

CD . 

43 ■ 



X 



L. Apanasevich, 1,2 J. Bacigalupi, 3 W. Baker, 4 M. Begel, 2 S. Blusk, 5 C. Bromberg, 1 P. Chang, 6 
B. Choudhary, 7 W. H. Chung, 5 L. de Barbara, 2 W. DeSoi, 2 W. Dhigosz, 6 J. Dunlea, 2 E. Engels, Jr., 5 
G. Fanourakis, 2 T. Ferbel, 2 J. Ftacnik, 2 D. Garelick, 6 G. Ginther, 2 M. Glaubman, 6 P. Gutierrez, 8 
K. Hartman, 9 J. Huston, 1 C. Johnstone, 4 V. Kapoor, 7 J. Kuehler, 8 C. Lirakis, 6 F. Lobkowicz, 2 ' * P. Lukens, 4 
J. Mansour, 2 A. Maul, 1 R. Miller, 1 B. Y. Oh, 9 G. Osborne, 2 D. Pellett, 3 E. Prebys, 2 R. Roser, 2 P. Shepard, 5 
R. Shivpuri, 7 D. Skow, 4 P. Slattery, 2 L. Sorrell, 1 D. Striley, 6 W. Toothacker, 9 * * S. M. Tripathi, 3 
N. Varelas, 2 D. Weerasundara, 5 J. J. Whitmore, 9 T. Yasuda, 6 C. Yosef, 1 M. Zielihski, 2 and V. Zutshi 7 

(Fcrmilab E706 Collaboration) 

1 Michigan State University, East Lansing, Michigan 48824 
2 University of Rochester, Rochester, New York 14627 
University of California-Davis, Davis, California 95616 
* Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
5 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
6 Northeastern University, Boston, Massachusetts 02115 
''University of Delhi, Delhi, India 110007 
8 University of Oklahoma, Norman, Oklahoma 73019 
9 Pennsylvania State University, University Park, Pennsylvania 16802 
(Dated: February 7, 2008) 

We present results on the production of direct photons, 7r°, and r\ mesons on nuclear targets 
at large transverse momenta {pt)- The data are from 530 and 800 GeV/c proton beams and 
515 GeV/c 7r~ beams incident upon copper and beryllium targets that span the kinematic range of 
1.0 < pt ^ 10 GeV/c at central rapidities. 

PACS numbers: 13.85.Qk, 12.38.Qk 



I. INTRODUCTION 

The study of inclusive particle production at large 
transverse momenta (pt) has yielded valuable infor- 
mation about perturbative quantum chromodynamics 
(PQCD), parton distribution functions (PDF), and frag- 
mentation functions of partons [1-5]. The use of nu- 
clear targets provides, in addition, information on parton 
and hadron rescattering and explores the time evolution 
of the collision. Since the discovery of the nuclear en- 
hancement of high-j»T single-particle production [6-8] , a 
large body of data has been accumulated to investigate 
nuclear-target effects in a wide variety of production pro- 
cesses, including those yielding single hadrons, dihadron 
pairs, Drell-Yan pairs, two-jet systems, and heavy fla- 
vors. Recent results from the RHIC program [9-12], in 
particular, have highlighted the differences between ini- 
tial and final-state effects in the nuclear environment. 
Many approaches have been developed to explain these 
data, which have included models for multiple-scattering, 
Fermi motion, modification of parton densities in the nu- 
clear medium, QCD higher-twist contributions, and new 
states of matter. 

We present the results of a high-statistics study of nu- 
clear effects in the inclusive production of direct photons, 
7r° and r\ mesons at large pt using data from Fermilab 
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experiment E706, and compare the results to predictions 
of a phenomenological model of nuclear effects [13]. 



II. APPARATUS 
A. Meson West spectrometer 

Fermilab E706 was a fixed-target experiment designed 
to measure the production of direct photons, neutral 
mesons, and associated particles at high-p^ [14-19]. The 
apparatus included a charged particle spectrometer and 
a large liquid argon calorimeter, as described below. Ad- 
ditional information about the Meson West spectrometer 
can be found in earlier papers [18, 20]. 

This paper reports on data from the two primary data 
runs of the experiment. During the 1990 run, the target 
consisted of two 0.8 mm thick copper foils followed by two 
pieces of beryllium (Fig. l:top). The upstream piece of 
beryllium was 3.7 cm long, while the length of the down- 
stream piece was 1.1 cm. In the 1991-1992 run, the target 
consisted of two 0.8 mm thick copper foils immediately 
upstream of a liquid hydrogen target [21], followed by 
a 2.54 cm long beryllium cylinder (Fig. Ebottom). The 
liquid hydrogen was contained in a 15.3 cm long mylar 
flask, which was supported in an evacuated volume with 
beryllium windows at each end (2.5 mm thickness up- 
stream and 2.8 mm thickness downstream). The target 
material is detailed in Table I. 

The charged particle spectrometer consisted of silicon 
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FIG. 1: Target configuration during the 1990 (top) and 1991- 
1992 (bottom) fixed target runs. 



TABLE I: Target fiducial lengths, densities, and number of 
interaction lengths by sample. The interaction lengths for 
the nuclear targets were obtained from Ref. [22]; for the H2 
target they were obtained from Ref. [23]. 



Target 


Beam 


Target Fiducial Measured Interaction 


Layout 


(GeV/c) 




Length 
(cm) 


Density 
(g/cm 3 ) 


Length 

(g/cm 2 ) 


1990 


515 71-- 


Be 


4.829 


1.848 


0.086 


Cu 


0.156 


8.810 


0.008 




515 n~ 








0.055 




530 p 


Be 


3.070 


1.855 


0.073 




800 p 








0.074 




515 7T~ 








0.008 


1991-1992 


530 p 
800 p 


Cu 


0.156 


8.810 


0.010 
0.010 




515 7T~ 








0.004 




530 p 


H 2 


14.500 


0.0705 


0.020 




800 p 








0.021 



microstrip detectors (SSDs) in the target region and mul- 
tiwire proportional chambers (PWCs) and straw tube 
drift chambers (STDCs) downstream of a large- aperture 
analysis magnet [18]. Six 3x3 cm 2 SSD planes were lo- 
cated upstream of the target region and used to recon- 
struct beam tracks. Two hybrid 5x5 cm 2 SSD planes 
(25 /jm pitch strips in the central 1 cm, 50 /im beyond) 
were located downstream of the target region. These 
were followed by eight 5x5 cm 2 SSD planes of 50 /jm 
pitch. The analysis dipole magnet imparted a 0.45 GeV/c 



Pt impulse in the horizontal plane to charged particles. 
Downstream track segments were measured by means of 
four stations of four views (XYUV) of 2.54 mm pitch 
PWCs and two stations of eight (4X4Y) layers of STDCs 
with tube diameters 1.03 cm (upstream station) and 
1.59 cm (downstream station) [24]. 

Photons were detected in a large, lead and liquid-argon 
sampling electromagnetic calorimeter (EMLAC), located 
9 m downstream of the target [20]. The EMLAC had a 
cylindrical geometry with an inner radius of 20 cm and 
an outer radius of 160 cm. The calorimeter had 33 lon- 
gitudinal cells read out in two sections: an 11 cell front 
section (8.5 radiation lengths) and a 22 cell back section 
(18 radiation lengths). Each longitudinal cell consisted 
of a 2 mm thick lead cathode (the first cathode was con- 
structed of aluminum), a double-sided copper-clad G-10 
radial (i?) anode board, a second 2 mm thick lead cath- 
ode, and a double-sided copper-clad G-10 azimuthal (<&) 
anode board. The 2.5 mm gaps between these layers were 
filled with liquid argon. The physical layout is illustrated 
in Fig. 2. 

The EMLAC readout was subdivided azimuthally into 
octants, each consisting of interleaved, finely segmented, 
radial and azimuthal views. This segmentation was real- 
ized by cutting the copper-cladding on the anode boards 
to form either radial or azimuthal strips. Signals from 
corresponding strips from all R (or $) anode boards 
in the front (or back) section of a given octant were 
jumpered together. The copper-cladding on the radial 
anode boards was cut into concentric strips centered on 
the nominal beam axis. The width of the strips on the 
first R board was 5.5 mm. The width of the strips on 
the following R boards increased slightly so that the ra- 
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dial geometry was projective relative to the target region. 
The azimuthal strips were split at a radius of 40 cm into 
inner and outer segments; each inner strip subtended an 
azimuthal angle of 7r/192 radians, while outer strips cov- 
ered 7r/384 radians. 

The spectrometer was located at the end of the Meson 
West beamline. The design of the beamline, primary tar- 
get, and primary beam dump were intended to minimize 
the rate of beam-halo muons incident upon the spectrom- 
eter. The beamline was capable of transporting either a 
primary (800 GeV/c) proton beam or unseparated sec- 
ondary particle beams of either polarity to the experi- 
mental hall. The beamline Cerenkov detector was used to 
identify the secondary beam particles [25]. This 43.4 m 
long helium-filled counter was located 100 m upstream 
of the experimental target. The positive secondary beam 
with mean momentum of 530 GeV/c was 97% protons. 
The negative secondary beam with mean momentum of 
515 GeV/c was 99% pions. 

At the end of the beamline was a 4.7 m long stack 
of steel surrounding the beam pipe and shadowing the 
EMLAC to absorb off-axis hadrons. A water tank was 
placed at the downstream end of this hadron shield to ab- 
sorb low-energy neutrons. Surrounding the hadron shield 
and neutron absorber were walls of scintillation counters 
(VW) to identify penetrating muons. There was one wall 
at the upstream end and two walls at the downstream 
end of the hadron absorber during the 1990 run. An 
additional wall was added to the upstream end of the 
hadron absorber prior to the 1991-1992 run. 



B. Trigger 

The E706 trigger selected interactions yielding high-pr 
showers in the EMLAC. The selection process involved 
several stages: beam and interaction definitions, a pre- 
trigger, and high-py trigger requirements [15, 18, 26]. A 
scintillator hodoscope, located 2 m upstream of the tar- 
get region, was used to detect beam particles, and reject 
interactions with more than one spatially isolated inci- 
dent particle. Additional scintillator with a 1 cm diame- 
ter central hole was located just downstream of the beam 
hodoscope, and served to reject interactions initiated by 
particles in the beam halo [27]. Two pairs of scintillator 
counters, mounted on the dipolc analysis magnet, were 
used to identify interactions in the target. To minimize 
potential confusion in the EMLAC due to out-of-time 
interactions, a filter was employed to reject interactions 
that occurred within 60 ns of one another. 

For those interactions that satisfied the beam and in- 
teraction requirements, the pt deposited in various re- 
gions of the EMLAC was evaluated by weighting the en- 
ergy signals from the EMLAC i?-channel amplifier fast 
outputs by a factor proportional to sin 6^, where 0i was 
the polar angle between the i th strip and the nominal 
beam axis. The PRETRIGGER hi requirement for a given 
octant was satisfied when the pt detected in either the in- 



ner 128 R channels or the outer R channels of that octant 
was greater than a threshold value. A pretrigger signal 
was issued only when there was no evidence in that oc- 
tant of substantial noise, significant px attributable to an 
earlier interaction, or incident beam-halo muon detected 
by the VW. 

Localized trigger groups were formed for each octant 
by clustering the i?-channel fast-outputs into 32 groups 
of 8 channels. Each adjacent pair of 8 channel groups 
formed a group-of-16 strips. If the pt detected in any 
of these groups-of-16 was above a specified high (or low) 
threshold, then a LOCAL hi (or LOCAL lo) signal was 
generated for that octant. A single local hi (or sin- 
gle local lo) trigger was generated if a LOCAL hi (or 
LOCAL lo) signal was generated in coincidence with the 
pretrigger HI in the same octant. 

Trigger decisions were also made based upon global 
energy depositions within an octant. A global lo sig- 
nal was generated if the total pt in an octant exceeded 
a threshold value. The LOCAL<g>GLOBAL LO trigger re- 
quired a coincidence of the pretrigger hi signal with 
GLOBAL LO and LOCAL LO signals from the same octant. 
The LOCAL LO requirement was included to suppress spu- 
rious global triggers due to coherent noise in the EMLAC. 

The single LOCAL LO and LOCAL<g>GLOBAL LO trig- 
gers were prescaled to keep them from dominating the 
trigger rate. Prescaled samples of beam, interaction, and 
pretrigger events were also recorded. 



III. ANALYSIS METHODS 

Data samples contributing to this analysis represent an 
integrated luminosity of 1.6 (6.8) pb -1 and 1.6 (6.5) pb _1 
for 530 and 800 GeV/c pCu (pBe) interactions, respec- 
tively, as well as 0.3 (1.4) pb" 1 for 515 GeV/c t^Cu 
(7r~Be) interactions. These samples were accumulated 
during the 1991-1992 run. Results reported in this pa- 
per also use 0.9 (6.1) pb _1 of 7r~Cu (7r~Be) data recorded 
during the 1990 run. The following subsections describe 
the analysis procedures and methods used to correct 
the data for losses from inefficiencies and selection bi- 
ases. Additional details can be found in our previous 
papers [15-18, 20]. 



A. Charged-particle reconstruction 

The two major aspects of the analysis procedure in- 
volved charged-particle and calorimeter-shower recon- 
struction (discussed in Sec. IIIB). The charged- track 
reconstruction algorithm produced track segments up- 
stream of the magnet using information from the SSDs, 
and downstream of the magnet using information from 
the PWCs and STDCs. These track segments were pro- 
jected to the center of the magnet and linked to form final 
tracks whose calculated charges and momenta were used 



for the physics analysis. The charged track reconstruc- 
tion is described in more detail elsewhere [18, 28]. The 
primary vertex reconstruction is described in Sec. HIE. 

B. Calorimeter shower reconstruction 

The readout of each EMLAC quadrant was divided 
into four regions: left and right R, and inner and outer $. 
Strip energies from clusters in each region were fit to the 
shape of an electromagnetic shower determined from de- 
tailed Monte Carlo simulations and isolated-shower data. 
These fits were used to evaluate the positions and ener- 
gies of the peaks in each region. Shower positions and 
energies were obtained by correlating peaks of approxi- 
mately the same energy in the R and $ regions within the 
same half octant. More complex algorithms were used to 
handle configurations with showers spanning multiple re- 
gions. The EMLAC readout was also subdivided longitu- 
dinally into front and back sections. This segmentation 
provided discrimination between showers generated by 
electromagnetically or hadronically interacting particles. 
Photons were defined as showers with at least 20% of the 
shower energy deposited in the front part of EMLAC, 
to reduce the backgrounds due to showers from hadronic 
interactions. Losses of photons due to this requirement 
were w 2%. A detailed event simulation was employed 
to correct for this and other effects including reconstruc- 
tion smearing and losses. An expanded discussion of the 
EMLAC reconstruction procedures and performance can 
be found elsewhere [20]. 
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FIG. 3: Comparison between data (histogram) and the de- 
tailed Monte Carlo simulation (•) for 515 GeV/c -k" beam on 
Cu for 77 combinations in the 7r° and r\ mass regions. 



in the 1990 run (1.08 in 1991-1992 run) and 1.19 per 
photon for the Cu target (1.16 in 1991-1992 run). 



D. Detector simulation 



C. Meson signals 

For this study, 7r° and r\ mesons were reconstructed via 
their 77 decay modes. Only those 77 combinations with 
energy asymmetry A 11 = \E 7l — E 72 \/(E 11 + E l2 ) < 0.75 
were considered to reduce uncertainties due to low en- 
ergy photons. The meson signals have been corrected for 
losses due to the energy asymmetry cut and the branch- 
ing fractions for the 77 decay modes [29] . 

Photons were required to be reconstructed within the 
fiducial region of the EMLAC to exclude areas with re- 
duced sensitivity. In addition, 77 combinations were re- 
stricted to the same octant to simplify the trigger anal- 
ysis. A simple ray-tracing Monte Carlo program was 
employed to determine the correction for these fiducial 
requirements. 

The correction for losses due to the conversion of pho- 
tons into e+e~ pairs was evaluated by projecting each 
reconstructed photon from the event vertex to the recon- 
structed position in the EMLAC. The radiation length of 
material traversed, up to the analysis magnet, was eval- 
uated based upon detailed detector descriptions. The 
photon conversion probability was evaluated and used to 
account for conversion losses. The average correction for 
conversion losses was 1.09 per photon for the Be target 



The Meson West spectrometer was modeled with a 
detailed GEANT [30] simulation (DGS). A preprocessor 
was used to convert GEANT information into the simu- 
lated hits and strip energies associated with the various 
detectors. The preprocessor simulated hardware effects, 
such as channel noise and gain variations. Monte Carlo 
events were then processed through the same reconstruc- 
tion software used for the analysis of the data. This 
technique accounted for inefficiencies and biases in the 
reconstruction algorithms. Reconstuction inefficiencies 
were relatively small over most of the kinematic range. 
More information on the detailed simulation of the Meson 
West spectrometer can be found elsewhere [15, 17, 31]. 

As inputs to the GEANT simulation, we employed sin- 
gle particle distributions, reconstructed data events, and 
HERWiG-generated [32] events. The HERWlG-generated 
7T°, 77, and direct-photon spectra were weighted in pt 
and rapidity to our measured results in an iterative fash- 
ion so that the final corrections were based on the data 
distributions rather than on the behavior of the physics 
generator. Figure 3 shows the 77 mass spectra in the 
7T° and 7] mass regions in comparison to the DGS re- 
sults for the 7T~Cu data at 515 GeV/c, and Fig. 4 shows 
an analogous plot for our higher statistics pBe data at 
530 GeV/c. In addition to providing evidence that the 
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FIG. 4: Comparison between data (histogram) and the de- 
tailed Monte Carlo simulation (•) for 530 GeV/c proton beam 
on Be for 77 combinations in the it and 7? mass regions. 



DGS simulated the EMLAC resolution well, the agree- 
ment between the levels of combinatorial background in- 
dicates that the DGS also provided a reasonable simu- 
lation of the underlying event structure. Since the DGS 
was tuned using our higher statistics Be data, Figs. 5-6 
also show the level of agreement achieved for this tar- 
get. Figure 5 shows a comparison between the DGS and 
the data for the sideband-subtracted energy asymmetry 
distribution for photons from 7r° decays. This figure illus- 
trates that the simulation accurately describes the losses 
of low-energy photons. 

A second Monte Carlo simulation of the detector 
(PMC) was used to cross check the detailed simulation 
and for studies that required large statistics. This simula- 
tion employed parametcrizations of physics cross sections 
and detector responses [17, 31, 33]. The inclusive 7r° and 
direct-photon cross sections were parameterized as two 
dimensional surfaces in px and rapidity [31]. The 77, u, 
and rj cross sections were parameterized using the mea- 
sured 77/71-° [15, 16], u/n [19], and rj'/r) [1, 34] ratios. 
Generated mesons were decayed into final state parti- 
cles; photons were smeared for energy and position res- 
olution [20]. A vertex was generated in the simulated 
target for every event. Photons were allowed to convert 
into e+e~ pairs; the energy of the resulting electrons was 
reduced using the GEANT function for bremsstrahlung ra- 
diation. Electron four-vectors were smeared for multiple 
scattering in the target and the resolution of the tracking 
system and adjusted for the magnet impulse. Figure 6 
displays a comparison between the PMC and the data in 
the 7T° and r\ mass regions and for the 7r° energy asym- 



FIG. 5: Comparison of A~ n , the energy asymmetry distri- 
bution for photons from ty° decays, in data (histogram) and 
the detailed Monte Carlo simulation (•) for the 800 GeV/c 
proton beam sample. Shown are the comparisons for two pr 
intervals, 4.0 < p T < 5.5 GeV/c and 5.5 <p T < 7.0 GeV/c. 



metry. The PMC provides an adequate characterization 
of the data. 



E. Vertex reconstruction 

The location of the interaction vertex was recon- 
structed using charged-particle tracks. Vertices were 
identified by means of an impact-parameter minimization 
technique [28]. A % 2 was defined for a given vertex po- 
sition using the impact parameters of the reconstructed 
tracks and their projection uncertainties. The vertex po- 
sition was found by minimizing this x 2 . Vertices were 
found in X and Y independently and correlated based on 
the difference in their positions along the nominal beam 
direction [Z axis). The Z position of the matched vertex 
was the weighted average of the Z positions found in X 
and Y . The reconstructed vertex positions are presented 
in Fig. 7 as functions of Z for the two target configura- 
tions. The beryllium, copper, and hydrogen targets are 
clearly visible, as are the SSDs and related support struc- 
tures. The average resolution for the Z location of the 
interaction vertex was « 300 /xm [18, 28]. 

The relative heights of the Cu and Be targets shown in 
Fig. 7 varies as a function of pr- This is clearly evident 
in Fig. 8, which compares two ir° samples: one acquired 
using the highly prescaled interaction trigger, and the 
other using the SINGLE LOCAL hi trigger. The inter- 
action triggered events are typically minimum-bias in 
character with low-px tt°'s. The number of primary ver- 
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FIG. 6: Comparison between data (•) and the parameterized 
Monte Carlo (histogram) from the 530 GeV/c proton beam 
sample: (top) 77 combinations in the 7r° and r\ mass regions; 
(bottom) asymmetry in energy for photons from 7r° mesons. 
The combinatorial background in the data has been removed 
from this comparison through a simple subtraction. The dis- 
tributions have been normalized to the same area. 



FIG. 8: Z positions of primary vertices for events containing 
7T° candidates with pr > 4.0 GeV/c acquired using the single 
LOCAL hi trigger (top) and for events containing 7r° candidates 
with pr > 0.5 GeV/c acquired using the interaction trigger 
(bottom). 
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FIG. 7: Z positions of primary vertices in the detailed Monte 
Carlo simulation (•) and the data (histogram) for events con- 
taining a 7T° candidate with pr > 4 GeV/c. 
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FIG. 9: Z positions of primary vertices in the vicinity of the 
Cu targets in the detailed Monte Carlo simulation (•) and the 
data (histogram), for events containing it candidates with 
pr > 4.0 GeV/c in the 800 GeV/c sample. Note the use of 
logarithmic scale. 
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tices scales as « A 2 / 3 , where A is the atomic weight of 
the target. The single local hi triggers are typically 
caused by hard-scatters that produce high-px ir°'s. The 
number of primary vertices in these events scale as « A 1 . 

The DGS was used for detailed studies of the ver- 
tex reconstruction. The transverse positions of ver- 
tices were chosen according to beam profiles observed in 
the data. Longitudinal positions were determined using 
Monte Carlo methods based upon the interaction lengths 
of the materials in the target region (Table I). DGS 
events were weighted to reproduce the relative number 
of vertices in the data. Results from the DGS compare 
favorably with the data in Fig. 7. This good agreement 
was particularly important for separating events with pri- 
mary vertices in copper from those in the upstream piece 
of beryllium. Figure 9 displays the longitudinal vertex 
distribution, focussing on this region in the 1991-1992 
target configuration. The shape of the tails in the data 
are well described by the DGS. 

The vertex reconstruction efficiency was evaluated us- 
ing the DGS [31]. Separate reconstruction efficiencies 
were evaluated for the Be, Cu, and H 2 targets. The re- 
construction probability was defined as the number of 
vertices reconstructed in each target's fiducial volume di- 
vided by the number of vertices generated in the fiducial 
volume. The reconstruction efficiency was the inverse of 
this probability. Defined in this manner, the reconstruc- 
tion efficiency also corrected for the longitudinal resolu- 
tion smearing of reconstructed vertices. 

Additional beam particles occasionally interacted in 
the target material within the data-capture timing win- 
dow of the tracking system. The extra tracks sometimes 
caused the vertex associated with the high-p-r interac- 
tion to be misidentified. The bias introduced by these 
rare events favored configurations where the low-pr in- 
teraction took place within the downstream piece of Be. 
This primarily affected interactions in the Cu and up- 
stream Be targets in the long 1991-1992 target config- 
uration because of the relatively poor vertex resolution 
in those targets compared to the downstream Be. This 
bias was investigated by comparing the ir° cross sections 
measured in 7r~Be interactions in the 1990 and the 1991- 
1992 runs, and by comparing ir° yields from the upstream 
and downstream Be pieces in the 1991-1992 target con- 
figuration. The number of Cu vertices were corrected for 
misidentifications arising from this source. The result- 
ing correction was 1.04 for the 1991-1992 7r~Cu sample 
at 515 GeV/c, 1.06 for the 530 GeV/c pCu sample, and 
1.12 for the 800 GeV/c pCu sample [35]. 

Each event in this analysis was required to have a re- 
constructed vertex in the target region. Longitudinal 
and transverse requirements were placed on vertices to 
define the data samples. The longitudinal cuts selected 
the target in which the incident beam interacted, while 
the transverse cuts ensured that the interaction occurred 
within the target material. 
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FIG. 10: A comparison of the direct-photon background as a 
function of pr as predicted by the parameterized Monte Carlo 
simulation (line) and the detailed Monte Carlo simulation (•). 
The background was normalized to the ir cross section. 



F. Direct photons 

The largest contribution to the direct-photon back- 
ground comes from electromagnetic decays of neutral 
hadrons, particularly 7r°'s and r/'s. For the purposes of 
the measurements reported here, a photon was a direct- 
photon candidate if it did not combine with another pho- 
ton in the same octant to form a 7r° with A 77 < 0.9 or 
an r\ with A 11 < 0.8. 

To suppress electrons, reconstructed showers were ex- 
cluded from the sample when charged-particle tracks 
pointed to within 1 cm of shower center. The correction 
for this criterion in the direct-photon analysis is « 1.01 
based upon studies of the impact of this requirement on 
reconstructed 7r°'s. 

The residual background from 7r°'s and rj's, as well as 
from other sources of background, was calculated using 
DGS samples that contained no generated direct photons 
(7t). A smooth fit of the -fb/iT in pt and rapidity was 
used to extract the direct-photon cross sections. The sys- 
tematic uncertainty in this background subtraction was 
estimated by varying the direct-photon definition as a 
function of the cut on ^4 77 . The direct-photon back- 
ground was also investigated using the PMC. Figure 10 
compares the direct-photon backgrounds estimated using 
the two Monte Carlo simulations. The close agreement 
provides additional confidence in our understanding of 
the background in our direct photon samples. Additional 
details are provided in Ref. [17]. 

Fits to 7h/7r° were only made for the beryllium target 
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FIG. 11: Difference between the copper and beryllium back- 
ground 7/71"° ratios versus pr for the 800 GeV/c proton beam. 
The dotted line represents a fit to the difference integrated 
over rapidity. 



due to relatively poor DGS statistics in the other tar- 
gets. However, 7b/7r° is expected to be slightly different 
for each target due to the different amounts of target ma- 
terial the photons must traverse. Therefore, a correction 
to 7f,/7r° was calculated using the PMC. The differences 
in 7b/7r° are shown for the 800 GeV/c proton beam data 
in Fig. 11. The other data samples have similar behav- 
ior. The target differences were fit as functions of px and 
rapidity for each incident beam and target configuration 
and applied as additive corrections to the nominal 76 / 7r° 
fit. 

The impact of the background can be determined by 
normalizing the direct-photon candidate spectrum from 
the data and the simulation to the measured 7r° cross 
section [15, 16]. This 7/71"° ratio is displayed in Fig. 12 
for all three incident beams for the copper target. The 
signal-to-background in all cases is large at high pt- 



G. Summary of systematic uncertainties 

The systematic uncertainties for the production of 
7r°'s, 37's, and direct photons measured in 7r~Cu andpCu 
interactions are similar to those detailed in Refs. [15-17]. 
The principal contributions to the systematic uncertainty 
arose from the following sources: normalization, calibra- 
tion of photon energy response and detector-resolution 
unsmearing, background subtraction, reconstruction effi- 
ciency, incident beam contamination (for the 515 GeV/c 
and 530 GeV/c secondary beams), beam halo muon re- 
jection, geometric acceptance, photon conversions, trig- 
ger response, and vertex finding. The additional vertex- 
finding uncertainty associated with the confusion induced 
by multiple beam particles interacting in the target was 
rs 2%. The total systematic uncertainties, combined in 
quadrature, are quoted with the cross sections in the ap- 
propriate tables. Note that some of these contributions 
to the systematic uncertainty (e.g. normalization) are 
strongly correlated between bins. 

Most of the experimental systematic uncertainties can- 
cel in the ratio of cross sections measured on different tar- 



pCu at 800 GeV/c 



♦ + 



! 7t"Cuat515GeV/c 



- • Data (unsubtracted) 
- Background 



5 7.5 10 5 7.5 10 5 7.5 10 

p T (GeV/c) 

FIG. 12: 7/71-° ratio for each of the three beam samples as 
a function of pr- Direct-photon candidates in the data are 
indicated by the points. The error bars represent statisti- 
cal contributions to the uncertainties. The line indicates the 
expected background and its width is representative of the 
associated systematic uncertainty. 



gets. The residual uncertainties are due to target-related 
systematics associated mainly with vertex identification. 
The total systematic uncertainty in the ratio of cross sec- 
tions measured on Cu to those on Be is ±3% for the 
800 GeV/c p beam sample, and ±2% for the 530 GeV/c 
p and 515 GeV/c 7T~ beam samples. The systematic un- 
certainty associated with the ratios of Be to H and Cu to 
H are ±4% for all samples after correcting for the effects 
of vertex misidentification as discussed in Sec. HIE. 



IV. RESULTS AND DISCUSSION 
A. Cross sections 

The invariant differential cross sections per nucleon 
for direct-photon, 7r°, and rj production from 530 and 
800 GeV/c p beams and 515 GeV/c 7T _ beam incident 
on copper are presented as functions of pt in Figs. 13 
through 15. Results from 530 GeV/c p and 515 GeV/c 
7r~ beams are averaged over the rapidity range —0.75 < 
Vcm < 0.75; results from the 800 GeV/c p beam are av- 
eraged over — 1.0 < y cm < 0.5. Data points are plotted 
at abscissa values that correspond to the average value 
of the cross section in each px bin, assuming local expo- 
nential px dependence [36] . The inclusive cross sections 
are tabulated in the Appendix. 
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FIG. 13: Invariant differential cross sections per nucleon for 
direct-photon, ir°, and r\ production as functions of pr, aver- 
aged over —0.75 < y cm < 0.75, from 530 GeV/c proton beam 
incident upon a copper target. 



FIG. 14: Invariant differential cross sections per nucleon for 
direct-photon, n°, and i] production as functions of pr, aver- 
aged over —1.0 < t/cm < 0.5, from 800 GeV/c proton beam 
incident upon a copper target. 



Ratios 



tion and agrees with the data within the uncertainties. 



E706 has previously reported results for direct-photon 
and 7T° production on beryllium and hydrogen tar- 
gets [15-17]. Since E706 is the only direct-photon exper- 
iment that used more than one nuclear target, our data 
provide a unique measurement of nuclear effects in direct- 
photon production. Figures 16 to 18 present the ratio 
of inclusive cross sections per nucleon measured on Cu 
to those measured on Be for n° mesons and direct pho- 
tons. (The Be/H ratios were presented in previous publi- 
cations [15-17], and are not reproduced here.) These ra- 
tios show clear evidence of nuclear enhancement in both 
7T° and direct-photon production and the n° data exhibit 
the decrease at high pr first noted by Cronin et al. [6, 7]. 
This behavior is generally attributed to the influence of 
multiple-parton scattering prior to the hard scatter [37]. 

Results from fits to constant ratios in restricted regions 
of pr have been overlaid on the data in Figs. 16 to 18, and 
summarized in Table II. Fits were made for 3.5 < pt < 
5.5 GeV/ c for the tt° 's and over the entire pr range for the 
direct photons. The difference in the ratio for tt° mesons 
and direct photons is significant and may be indicative 
of differing roles of initial and final-state effects, since 
direct photons are not expected to be strongly impacted 
by final-state nuclear effects. 

Expectations from a theoretical prediction for nuclear 
enhancement in the 515 GeV/c ir~ direct-photon sam- 
ple [38] have been overlaid on Fig. 18. This calculation 
predicts a slight enhancement in direct-photon produc- 



C. Comparisons with hijing 

The Cu to Be cross-section ratios are compared with 
results from the hijing Monte Carlo event generator in 
Figs. 19 and 20. hijing is a program designed to simulate 
particle production in pp, pA, and AA collisions [13]. It 
was necessary to normalize the HIJING results in Figs. 19 
and 20 to the data. However, the shapes of the curves, in 
the case of the proton beams, are in good agreement with 
the data. Renewed interest in the amount of nuclear en- 
hancement as a function of rapidity has been generated 
by recent BRAHMS measurements [39] and correspond- 
ing results from the other RHIC experiments [9-11]. Our 
high-statistics cross-section measurements in the central 
rapidity region can be used to tune theoretical models 
developed to describe the RHIC environment. The ra- 
pidity dependence of the Cu to Be cross-section ratios 



TABLE II: Results of fits of the Cu to Be ratios for tt° and 
direct-photon production. 



Sample 


7T° 

3.5 <p T < 5.5 GeV/c 


direct photon 
3.5 <p T < 8.0 GeV/c 


530 GeV/c p 
800 GeV/c p 
515 GcV/c 7T- 


1.271 ±0.016 ±0.025 
1.283 ± 0.025 ±0.038 
1.237 ±0.015 ±0.025 


1.103 ±0.032 ±0.022 
1.043 ±0.032 ±0.031 
1.083 ±0.024 ±0.022 
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FIG. 15: Invariant differential cross sections per nucleon for 
direct-photon, tv°, and r\ production as functions of pr, av- 
eraged over —0.75 < y cm < 0.75, for 515 GeV/c tt~ beam 
incident upon a copper target. 
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FIG. 16: The ratio of inclusive it and direct-photon produc- 
tion cross sections per nucleon in pCu to those in pBe colli- 
sions at 530 GeV/c. Simple straight line fits to regions with 
relatively flat distributions have been overlaid on the data. 
The error bars represent only statistical contributions to the 
uncertainties. Systematic uncertainties are indicated by the 
shaded region associated with the fit. 



FIG. 17: The ratio of inclusive it and direct-photon produc- 
tion cross sections per nucleon in pCu to those in pBe colli- 
sions at 800 GeV/c. Simple straight line fits to regions with 
relatively flat distributions have been overlaid on the data. 
The error bars represent only statistical contributions to the 
uncertainties. Systematic uncertainties are indicated by the 
shaded region associated with the fit. 



are compared to expectations from hijing in Figs. 21 to 
23. hijing does not describe the rapidity dependence of 
the 7T° data for the incident proton beams; the hijing 
results are generally peaked towards backward rapidities 
(like the BRAHMS data), whereas our data are relatively 
independent of rapidity, hijing provides a better descrip- 
tion of our direct-photon data. 



V. CONCLUSIONS 

We have measured the invariant differential cross sec- 
tion per nucleon for direct-photon, 7r°, and rj production 
from 515 GeV/c ir~ beam and 800 and 530 GeV/c pro- 
ton beams incident on copper as a function of px and 
y cm _. These data span the kinematic range 1.0 < pt 
10 GeV/c and central rapidities. 

Ratios of these production cross sections to our previ- 
ously published measurements on a beryllium target [15- 
17] show a strong nuclear enhancement for 7r° mesons 
and a smaller, but significant, enhancement for direct 
photons. We compare these measurements with expec- 
tations from a theoretical calculation (for direct-photon 
production for incident 7r~ beam) and with the results of 
a Monte Carlo event generator, hijing, dedicated to the 
simulation of nuclear effects, hijing yields a good de- 
scription of the shape of the px dependence of the Cu to 
Be ratios for both 7r°'s and direct photons in the incident 
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FIG. 18: The ratio of inclusive it and direct-photon produc- 
tion cross sections per nucleon in 7r~Cu to those in 7r~Be 
collisions at 515 GeV/c. A theoretical prediction for direct- 
photon production from X. Guo and J. Qiu [38] is overlaid on 
the data (dotted curve). Simple straight line fits to regions 
with relatively fiat distributions have also been overlaid on the 
data. The error bars represent only statistical contributions 
to the uncertainties. Systematic uncertainties are indicated 
by the shaded region associated with the fit. 
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FIG. 19: The ratio of inclusive 7r° production cross sections 
per nucleon on Cu to those on Be, compared with predictions 
from hijing scaled to the data. The error bars represent only 
statistical contributions to the uncertainties and v represents 
the number of degrees of freedom of the fit. 



FIG. 20: The ratio of inclusive direct-photon production cross 
sections per nucleon on Cu to those on Be, compared with 
predictions from hijing scaled to the data. The error bars 
represent only statistical contributions to the uncertainties 
and v represents the number of degrees of freedom of the fit. 
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FIG. 21: The rapidity dependence of the ratio of inclusive 
cross sections per nucleon on Cu to those on Be for direct- 
photon and 7T° production in the 530 GeV/c p beam, com- 
pared with predictions from hijing scaled to the data. The 
error bars represent only statistical contributions to the un- 
certainties and v represents the number of degrees of freedom 
of the fit. 
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proton beam data samples, hijing also describes the ra- 
pidity dependence of direct photon production in those 
samples. However, hijing provides a relatively poor de- 
scription of the rapidity dependence of tt° meson produc- 
tion. 



FIG. 22: The rapidity dependence of the ratio of inclusive 
cross sections per nucleon on Cu to those on Be for direct- 
photon and 7T° production in the 800 GeV/c p beam, com- 
pared with predictions from hijing scaled to the data. The 
error bars represent only statistical contributions to the un- 
certainties and v represents the number of degrees of freedom 
of the fit. 
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FIG. 23: The rapidity dependence of the ratio of inclusive 
cross sections per nucleon on Cu to those on Be for direct- 
photon and 7r production in the 515 GeV/c ir~~ beam, com- 
pared with predictions from hijing scaled to the data. The 
error bars represent only statistical contributions to the un- 
certainties and v represents the number of degrees of freedom 
of the fit. 



Acknowledgments 



We thank the U. S. Department of Energy, the Na- 
tional Science Foundation, including its Office of Inter- 
national Programs, and the Universities Grants Commis- 
sion of India, for their support of this research. The staff 
and management of Fermilab are thanked for their ef- 
forts in making available the beam and computing fa- 
cilities that made this work possible. We are pleased 
to acknowledge the contributions of our colleagues on 
Fermilab experiment E672. We acknowledge the contri- 
butions of the following colleagues to the operation and 
upgrade of the Meson West spectrometer: W. Dickerson 
and E. Pothier from Northeastern University; J. T. An- 
derson, E. Barsotti Jr., H. Koecher, P. Madsen, D. Pc- 
travick, R. Tokarek, J. Tweed, D. Allspach, J. Urbin, and 
the cryo crews from Fermi National Accelerator Labora- 
tory; T. Haelen, C. Benson, L. Kuntz, and D. Ruggiero 
from the University of Rochester; the technical staffs of 
Michigan State University and Pennsylvania State Uni- 
versity for the construction of the straw tubes and of the 
University of Pittsburgh for the silicon detectors. We 
thank the following commissioning run collaborators for 
their invaluable contributions to the hardware and soft- 
ware infrastructure of the original Meson West spectrom- 
eter: G. Alverson, G. Ballocchi, R. Benson, D. Berg, 
D. Brown, D. Carey, T. Chand, C. Chandlce, S. Easo, 
W. Faissler, G. Glass, I. Kourbanis, A. Lanaro, C. Nel- 
son Jr., D. Orris, B. Rajaram, K. Ruddick, A. Sinanidis, 
and G. Wu. We also thank X.-N. Wang for several helpful 
discussions. 



13 



[1] W. M. Geist et al., Phys. Rept. 197, 263 (1990). 
[2] N. A. McCubbin, Rep. Prog. Phys. 44, 65 (1981). 
[3] J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 
[4] T. Ferbel and W. R. Molzon, Rev. Mod. Phys. 56, 181 
(1984). 

[5] L. Apanasevich et al., Phys. Rev. D59, 074007 (1999). 
[6] J. W. Cronin et al., Phys. Rev. Dll, 3105 (1975). 
[7] D. Antreasyan et al., Phys. Rev. D19, 764 (1979). 
[8] H. Frisch et al., Phys. Rev. D27, 1001 (1983). 
[9] J. Adams et al. (STAR), nucl-ex/0501009. 
[10] K. Adcox et al. (PHENIX), nucl-ex/0410003. 
[11] B. B. Back et al. (PHOBOS), nucl-ex/0410022. 
[12] I. Arsene et al. (BRAHMS), nucl-ex/0410020. 
[13] X.-N. Wang, Nucl. Phys. A661, 609 (1999); Private com- 
munication with Xin-Nian Wang regarding the use of this 
program at these beam energies. 
[14] L. Apanasevich et al. (E706), Phys. Rev. Lett. 81, 2642 
(1998). 

[15] L. Apanasevich et al. (E706), Phys. Rev. D68, 052001 

(2003) . A computational error in the correction of mea- 
sured cross sections for contributions from minority par- 
ticles in secondary beams affected the reported cross sec- 
tions for 7T° and r\ production by 530 GeV/c protons. 
As a result, the cross sections listed in the tables for 
530 GeV/c protons in the referenced publication should 
be increased by 3.4%, a change that falls well within the 
quoted systematic uncertainties. 

[16] L. Apanasevich et al. (E706), Phys. Rev. D69, 032003 

(2004) . 

[17] L. Apanasevich et al. (E706), Phys. Rev. D70, 092009 

(2004) . 

[18] L. Apanasevich et al. (E706), Phys. Rev. D56, 1391 
(1997). 

[19] L. Apanasevich et al. (E706), hep-ex/0004012. 
[20] L. Apanasevich et al. (E706), Nucl. Inst. & Meth. A417, 
50 (1998). 

[21] D. Allspach et al., Adv. Cryog. Eng. 37, 1495 (1991). 

[22] A. S. Carroll et al., Phys. Lett. 80B, 319 (1979). 

[23] A. Baldini et al., in Total Cross-Sections for Reactions of 

High Energy Particles, edited by H. Schopper (Springer- 

Verlag, New York, 1988), vol. 12. 
[24] C. Bromberg et al. (E706), Nucl. Inst. & Meth. A307, 

292 (1991). 

[25] D. Striley, Ph.D. thesis, University of Missouri at 
Columbia (1996). 

[26] L. Sorrell, The E706 Trigger System, E706 Note 201, 
1994 (unpublished); Ph.D. thesis, Michigan State Uni- 
versity, 1995; G. Osborne, Ph.D. thesis, University of 
Rochester, 1996. 

[27] During the 1990 run, the beam hole definition was im- 
plemented using a single scintillation counter. An array 
of four scintillation counters was used for the beam hole 
definition during the 1991-92 run. 

[28] S. Blusk, Ph.D. thesis, University of Pittsburgh (1995). 

[29] D. E. Groom et al. (Particle Data Group), Eur. Phys. J. 
C15, 1 (2000). 

[30] F. Carminati and others, GEANT: Detector Descrip- 
tion and Simulation Tool, CERN Program Library Long 
Writeup W5013, 1993. 

[31] L. Apanasevich, Ph.D. thesis, Michigan State University 

(2005) . 



[32] G. Marchesini et al., Comput. Phys. Commun. 67, 465 
(1992), herwig v5.6. 

[33] M. Begel, Ph.D. thesis, University of Rochester (1999). 

[34] M. Diakonou et al. (R806), Phys. Lett. 89B, 432 (1980). 

[35] Our previously published cross sections on hydrogen were 
not corrected for these misidentified vertices (the corre- 
sponding corrections in the case of Be were negligible). 
The applicable correction factors are 1.01 for the 1991- 
92 tt~p sample at 515 GeV/c, 1.02 for the 530 GeV/c pp 
sample, and 1.06 for the 800 GeV/c pp sample. In each 
case, the additional correction falls within our published 
systematic uncertainty. 

[36] G. D. Lafferty and T. R. Wyatt, Nucl. Instrum. Methods 
Phys. Res. A355, 541 (1995). 

[37] X.-N. Wang, Phys. Rept. 280, 287 (1997). 

[38] X. Guo and J. Qiu, Phys. Rev. D53, 6144 (1996). 

[39] I. Arsene et al. (BRAHMS), Phys. Rev. Lett. 93, 242303 
(2004). 



14 



APPENDIX: TABULATED CROSS SECTIONS 

In this appendix, we present tables of the measured invariant differential cross sections for direct photon, ir° , and 
r\ production on Cu targets as functions of pt- In these tables, the first uncertainty is statistical and the second is 
systematic. In the case of the lowest two pr bins for the n° measurement, the statistical and systematic uncertainties 
have been combined because of the large correlation between them. 

TABLE III: Invariant differential cross sections (Eda/d 3 p) per nucleon for direct-photon production in pCu collisions at 800 
and 530 GeV/c, and 7r~Cu collisions at 515 GeV/c as functions of pr- 



PT 


pCu at 530 GeV/c 


pCu at 800 GeV/c 


tt"Cu at 515 GeV/c 






-0.75 < 2/cm < 0.75 


-1.0 < ?/cm < 0.5 


-0.75 < y C m < 0.75 


(GeV/c) 


[nb/(GeV/c) 2 ] 


[nb/(GeV/c) 2 ] 


[nb/(GcV/c) 2 ] 


3.50 - 


3.75 


3.11 ± 0.45 ±0.63 


3.01 ± 0.59 ±0.63 


1.84 ± 0.28 ±0.34 


3.75 - 


4.00 


1.21 ± 0.15 ±0.22 


1.87 ± 0.23 ±0.36 


1.065 ± 0.095 ±0.18 


4.00 - 


4.25 


0.583 ± 0.048 ± 0.099 


0.989 ± 0.089 ±0.17 


0.585 ± 0.022 ± 0.093 


4.25 - 


4.50 


0.282 ± 0.011 ± 0.045 


0.523 ± 0.026 ±0.086 


0.361 ± 0.015 ±0.053 


4.50 - 


4.75 


0.1516 ± 0.0074 ±0.023 


0.299 ± 0.017 ±0.046 


0.1783 ± 0.0092 ±0.025 


4.75 - 


5.00 


0.0865 ± 0.0052 ±0.012 


0.182 ± 0.011 ±0.026 


0.1124 ± 0.0066 ± 0.015 






[ P b/(Gcv/ c n 


[pb/(GeV/ C n 


[pb/(GeV/cr] 


5.00 - 


5.25 


50.2 ± 3.5 ±6.7 


110.6 ± 7.6 ±15 


81.0 ± 4.8 ± 10 


5.25 - 


5.50 


27.7 ± 2.5 ± 3.6 


59.9 ± 5.5 ±7.9 


40.4 ± 3.4 ±4.9 


5.50 - 


5.75 


15.6 ± 1.8 ±2.0 


37.0 ± 4.1 ±4.7 


32.3 ± 2.7 ±3.9 


5.75 - 


6.00 


11.2 ± 1.3 ± 1.4 


21.4 ± 2.9 ±2.7 


23.1 ± 2.1 ± 2.7 


6.00 - 


6.50 


4.99 ± 0.62 ±0.61 


15.4 ± 1.3 ±1.9 


9.74 ± 0.95 ±1.1 


6.50 - 


7.00 


1.93 ± 0.34 ±0.24 


5.43 ± 0.72 ±0.64 


5.30 ± 0.60 ±0.61 


7.00 - 


8.00 


0.278 ± 0.094 ± 0.035 


2.01 ± 0.27 ±0.23 


1.84 ± 0.23 ±0.22 


8.00 - 


10.00 


0.086 ± 0.029 ±0.012 


0.136 ± 0.055 ±0.016 


0.197 ± 0.052 ±0.024 


10.00 - 


12.00 




0.0077 ± 0.0077 ± 0.0009 


0.0124 ± 0.0097 ± 0.0016 
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TABLE IV: Invariant differential cross sections (Eda/d 3 p) per nucleon for n° production in pCu collisions at 800 and 530 GeV/c, 
and 7r~Cu collisions at 515 GeV/c as functions of pr- 





Pt 


pCu at 530 GeV/c 


pCu at 800 GeV/c 


tt-Cu at 515 GeV/c 






-0.75 < y C m < 0.75 


-1.0 < I/cm < 0.5 


-0.75 < y C m < 0.75 


(GeV/c) 


[/xb/(GeV/c) 2 ] 


[/xb/(GeV/c) 2 ] 


[ M b/(GeV/c) 2 ] 


l.UU 


1 /in 
— 1.4U 


315 ± 50 


293 ± 67 


166 ± 32 


1.4U 


i on 
- 1.8U 


43.2 ± 9.9 


79 ± 16 


23.9 ± 6.2 


i on 
l.oU 


o on 


7.3 ± 2.6 ±0.8 


18.5 ± 4.1 ±2.0 


8.72 ± 0.88 ± 1.0 


o on 


o a n 
— Z.4U 


2.64 ± 0.19 ±0.28 


3.77 ± 0.40 ±0.42 


2.015 ± 0.094 ±0.24 


o /in 
z.4U 


o fin 
— z.DU 


1.38 ± 0.23 ±0.15 


2.00 ± 0.20 ±0.22 


0.973 ± 0.057 ±0.11 


o fin 


o on 
— z.oU 


0.663 ± 0.041 ±0.070 


0.797 ± 0.053 ±0.088 


0.505 ± 0.023 ±0.058 


2.80 


- 3.00 


0.2723 ± 0.0093 ± 0.029 


0.476 ± 0.014 ± 0.052 


0.244 ± 0.012 ±0.028 


3.00 


- 3.20 


0.1453 ± 0.0048 ±0.016 


0.2287 ± 0.0082 ±0.025 


0.1417 ± 0.0062 ±0.016 






[nb/(GeV/c) 2 l 

L / V / / J 


[nb/(GeV/c) 2 l 

L / V / / J 


fnb/(GeV/c) 2 l 


q on 
o.zU 


q /in 
— o.4U 


77.5 ± 3.2 ±8.2 


132.5 ± 5.3 ± 15 


69.2 ± 2.7 ±7.7 


q /in 


Q fin 
— o.DU 


42.0 ± 2.1 ±4.5 


70.1 ± 3.5 ±7.7 


39.1 ± 2.0 ±4.3 


q fin 
O.OU 


q sn 
— o.oU 


20.0 ± 1.2 ± 2.1 


41.2 ± 2.6 ±4.5 


21.20 ± 0.92 ± 2.4 


O.oU 


/i nn 

— 4.UU 


11.53 ± 0.34 ± 1.2 


23.1 ± 1.4 ±2.5 


12.27 ± 0.37 ± 1.4 


/i nn 
4.UU 


/i on 

— 4.ZU 


6.411 ± 0.067 ±0.69 


12.74 ± 0.68 ±1.4 


6.947 ± 0.075 ±0.77 


/i on 
4.zU 


/i /in 
— 4.4U 


3.579 ± 0.043 ± 0.39 


6.89 ± 0.15 ±0.76 


4.021 ± 0.054 ±0.44 


A Af\ 
4.4U 


a fin 
— 4.DU 


2.002 ± 0.030 ±0.22 


4.15 ± 0.12 ±0.46 


2.389 ± 0.035 ±0.26 


A fin 
4.0U 


/i en 
— 4.oU 


1.146 ± 0.020 ± 0.13 


2.528 ± 0.075 ± 0.28 


1.438 ± 0.024 ±0.16 


a sn 

4.oU 


^ nn 
- o.uu 


0.639 ± 0.014 ±0.070 


1.426 ± 0.041 ±0.16 


0.892 ± 0.018 ±0.098 






[pb/(GeV/cfl 


fpb/(GeV/c) 2 l 


[pb/(GcV/c) 2 l 


5.00 


- 5.20 


392 ± 11 ±43 


931 ± 31 ± 110 


514 ± 13 ± 57 


5.20 


- 5.40 


215.8 ± 7.6 ±24 


582 ± 23 ± 66 


320 ± 11 ± 35 


5.40 


- 5.60 


134.3 ± 5.9 ± 15 


374 ± 18 ± 42 


219.2 ± 8.3 ± 24 


5.60 


- 5.80 


72.2 ± 4.3 ± 8.1 


258 ± 14 ± 29 


130.9 ± 6.7 ± 15 


5.80 


- 6.00 


44.4 ± 3.2 ± 5.0 


143.6 ± 9.3 ± 16 


79.7± 4.7 ±8.9 


6.00 


- 6.25 


25.2 ± 2.0 ± 2.9 


88.7 ± 5.2 ± 10 


51.2 ± 3.3 ±5.8 


6.25 


- 6.50 


16.6± 1.7 ± 1.9 


53.9 ± 3.7 ±6.3 


27.8 ± 2.4 ±3.2 


6.50 


- 6.75 


6.3 ± 1.0 ±0.7 


33.8 ± 3.1 ±3.9 


15. 3± 1.7 ± 1.7 


6.75 


- 7.00 


3.71 ± 0.71 ±0.44 


19.1 ± 1.9 ±2.2 


9.5± 1.4 ± 1.1 


7.00 


- 7.50 


1.96 ± 0.38 ±0.23 


9.5± 1.0 ± 1.1 


4.39 ± 0.63 ±0.51 


7.50 


- 8.00 


0.54 ± 0.18 ±0.07 


2.10 ± 0.46 ±0.25 


1.73 ± 0.44 ±0.21 


8.00 


- 9.00 


0.149 ± 0.061 ±0.019 


0.95 ± 0.26 ±0.12 


0.43 ± 0.17 ±0.05 


9.00 


- 10.00 






0.067 ± 0.047 ±0.009 



TABLE V: Invariant differential cross sections (Eda/d 3 p) per nucleon for r\ production in pCu collisions at 800 and 530 GeV/c, 
and 7r~Cu collisions at 515 GeV/c as functions of pr- 





Pt 


pCu at 530 GeV/c 


pCu at 800 GeV/c 


tt-Cu at 515 GeV/c 






-0.75 < ?/cm < 0.75 


-1.0 < J/cm < 0.5 


-0.75 < ?/cm < 0.75 


(GeV/c) 


[nb/(GeV/c) 2 ] 


[nb/(GeV/c) 2 ] 


[nb/(GeV/c) 2 ] 


3.00 


- 3.50 


53 ± 12 ± 6.0 


82 ± 25 ± 10 


44 ± 12 ± 6.0 


3.50 


- 4.00 


7.9 ± 2.2 ± 0.9 


9.9 ± 5.1 ± 1.2 


10.6 ± 1.6 ± 1.3 


4.00 


- 4.50 


1.68 ± 0.17 ±0.19 


4.78 ± 0.47 ±0.56 


2.21 ± 0.17 ±0.26 


4.50 


- 5.00 


0.434 ± 0.043 ± 0.050 


1.06 ± 0.12 ±0.12 


0.730 ± 0.049 ± 0.084 






[pb/(GeV/c)1 


[pb/(GeV/c)1 


[ P b/(GeV/c)1 


5.00 


- 5.50 


120 ± 13 ± 14 


274 ± 37 ± 32 


203 ± 17 ± 24 


5.50 


- 6.00 


41.0 ± 6.6 ±4.9 


123 ± 23 ± 15 


56.6 ± 8.0 ±6.6 


6.00 


- 7.00 


7.5 ± 1.4 ±0.9 


15.8 ± 4.7 ± 1.9 


18.2 ± 2.3 ±2.2 


7.00 


- 8.00 


1.01 ± 0.35 ±0.13 


3.0 ± 1.5 ±0.4 


2.66 ± 0.76 ± 0.33 


8.00 


- 10.00 0.087 ± 0.084 ±0.012 


0.15 ± 0.41 ±0.02 


0.29 ± 0.16 ±0.04 



